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Abstract 

A matrix norm gives an upper bound on the spectral radius of a matrix. Knowledge on the 
location of the dominant eigenvector also leads to upper bound of the second eigenvalue. We show 
how this technique can be used to prove that certain semi-group of matrices arising from continued 
fractions have a Pisot spectrum: namely for all matrices in this semi-group all eigenvalues except 
the dominant one is smaller than one in absolute value. 


1 Introduction 


A dominant eigenvalue of a real square matrix is an eigenvalue of maximum modulus. We call a 
square matrix Pisot if it has non-negative integer entries, its dominant eigenvalue is simple and 
all eigenvalues different from the dominant one have absolute values less than one. We prove that 
several monoids of non-negative matrices enjoy the property of all being Pisot. 

Our first family of matrices is related to the so called fully subtractive (multidimensional) con¬ 


tinued fraction algorithm. For an integer d > 2 we define for each k = 1, 
by 

fAk) \ _ / 1 a j = k or i = j, 

\^FS,d)A — 


, d the matrix A 


(fc) 

FS,d 


0 otherwise 


For d = 3 this boils down to the three matrices 


4 ( 1 ) _ 

^FS,3 ~ 



4(2) 

^FS,3 



4 ( 3 ) _ 

^FS,3 ~ 



All non-degenerate products of the matrices satisfy the Pisot property. 


Theorem 1. Let A = ^ be a product of the fully subtractive matrices in dimen¬ 

sion d. Then the matrix A is primitive if and only if all letters {1,..., d} appear in the sequence 
(fi, 12 ,..., in)- Moreover, if the matrix A is primitive then it is Pisot. 

Recall that a non-negative square matrix A is primitive if there exists a positive integer n so 
that A" has all its entries positive. The case d = 3 of Theorem]^ was proved in |ArIt01| . The 
authors used an induction on characteristic polynomials and our approach is radically different. 

The same result holds for another set of 3 x 3 matrices related to the Brun multidimensional 
continued fractions. Let 





Theorem 2. Let A = A^^)} A^^J ... A^"^ be a product of the 3x3 Brun matrices. Then, B is 
primitive if and only if the matrix B'A) appears in the product. Moreover, if B is primitive then it 
is Pisot. 

This result was already known since the work of Brun |Br57| . 

The proofs of Theorem and are elementary and uses the following inequality. Given a non¬ 
negative primitive d x d matrix A and its Perron-Frobenius eigenvector v G 
A 2 of its second largest eigenvalue satisfies 

^ IIAxll 

A 2 < sup -n—ip- 


the absolute value 


1 







where ||.|| is any norm on In our proof, the information we have on the localization of the 
Perron-Frobenius eigenvector comes from the dynamical systems induced by the matrices; in other 
words the fully subtractive and Brun continued fraction algorithms. 

From a diophantine approximation point of view, the Pisot property is particularly interesting 
because it provides the so called exponential convergence of the continued fraction expansion for 
almost every vectors (see |LaM] l. We show that the above results naturally extends to this situation 
in Section [6l 

Beyond continued fractions, Pisot matrices are of special interest in substitutive dynamical 
systems. More precisely, replacing matrices with so called substitutions, the associated dynamical 
systems admit d — 1 eigenvalues and in many cases the dynamical system can be proved to have 
purely discrete spectrum (see |Py02 Chapter 7). This was the main motivation for the study of 
the fully subtractive matrices in |ArIt01| . 

2 Fully subtractive and Brun continued fractions 

Let us consider a finite or countable set A that we call alphabet and for each i G A a matrix 
G SL{d,Z). We already saw two examples of this with the fully subtractive algorithm where 
■^FS,d = {1) 2,..., d} and for Brun algorithm where Asr = {1, 2, 3}. 

To the data (A, we associate the set of infinite words A = the shift map T : A —>■ 

A and a cocycle 

Vcc G A, Vn > 0, A„{x) = 

The maps An : A —>■ SL(d, Z) satisfy the so called cocycle property. Am+ni.x') — Ani{x)An{T x). 

Definition 3. Let be a set of matrices in SL(d,Z) where A is a finite or countable 

alphabet. We say that a set D C P(M‘^) is adapted to these matrices if it is non-empty, it is the 
closure of its interior and for all i G A we have A^'^'^D C D. 

For example PP{M.'^) is always adapted. But we will be interested in the somewhat smallest 
adapted set in order to localize the dominant eigenvector. 

Given and D C P(K‘^) adapted, we define = A^'^^D and more generally for a 

finite word w = iAi ■ ■ ■ in-i we define ... A^'‘'^-^'>D. Note that for any w the set 

is not empty. Given an infinite word x = XqXi ... G A we also set Dn{x) = An{x)D and 

Dooix) = nn>0 ^ri(.x). 

Let A G SL(d, Z) and let A be its spectral radius. Consider its Jordan decomposition over C and 
the Jordan blocks associated with an eigenvalue of modulus A and being of maximal dimension. 
To each of these maximal Jordan block is associated exactly one eigenvector Vi. The dominant 
eigenspace of Al is n (Cui 0 Cv 2 © ... Cu^). We have the following elementary result. 

Lemma 4. Let be a finite or countable set of matrices in SL((i, Z) and let D be adapted. 

Let X = xqXi ... G a. be an infinite word over A. Then for any n, the set Dn{x) contains a basis 
of the dominant eigenspace of An{x). 

We omit the proof that only uses the fact that the maximum growth of ||A"u|j is A"n^ where 
k is the maximal dimension of a Jordan block associated with an eigenvalue of maximal modulus 
of A. 

Let 

Dps.d = {(a^i,. ..,Xd)G P(K+) : Vi, j, fc Xi < xj 0 Xk} 

and 

L>Br = {ix,y,z) G P(M+) : x> y > z}. 

Then it is easily seen that Dps.d is adapted for the fully subtractive matrices in dimension d and 
Dbt is adapted for the Brun matrices. In figures and one can see the projective picture of the 
domains ZJV), and Note that in these cases, the domains are disjoint but that it is 

not a requirement in our definition. Moreover, one can see that in the Brun case the form a 
partition while it is not the case for the fully subtractive. 

If the D(®) are disjoint one can define a continued fraction algorithm as follows. One defines a 
partial map / : D D by setting fx = (A(®))“^x on 
One can compute that for Brun one has 

fBr{x) = sort(x - y,y,z) 
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.( 1 ) _ 
^FS,3 — 



FS,3 


.0 1 


4(3) _ 

^FS,3 ~ 



( 110 ) 


Figure 1: Fully subtractive partition of the domains with d = 3. 

(Ill) 

4(2) _ 

^Br ~ 




yig) = I 1 0 0 
.0 1 0 > 


'1 0 1^ 


( 100 ) 



( 110 ) 


Figure 2: The matrices and domains for the Brun algorithm. 


where sort : P(Mi)_) D is the map which permutes the coordinates in order to sort them. While 
for the fully subtractive one has 


fpsA^) = X2- Xi,..., Xi-i - Xi, Xi, Xi+i -Xi,...,Xd-Xi) if Xi = min(a;i, ...,Xd). 


3 Strategy 

The proofs of Theoremsandfollow a general strategy that we describe now. We let ||.|j be the 
L°° norm on and the associated operator norm on matrices. That is for a vector v and a matrix 
A 

||u|| = max(|ui|,..., |ud|) and \\A\\ = max |a„|. 

In this section the norm used on has no importance. But it turns out that, to apply the results 
to continued fraction algorithms, the most convenient one was always the norm. 

To a non-zero vector v in we associate its dual hyperplane Hy = {z S {v, z) = 0}. Given 

a non zero vector v in we define the following semi-norm on d x d matrices 

\\B\\v= sup = max ||Sz||. 


More generally, if A C is a cone, we define 


||B||a= sup ||S||„ 
«eP(A) 


Let (A(*))jg^ be a finite or countable set of matrices as in Section]^ Let also A = and 
D C P(K‘^) be adapted. Recall that Dn{x) = An{x)D = particular, Dq{x) = D 

and Di{x) = and that D^{x) = 


Lemma 5. Let (A(®))ig^ be a finite or countable set of matrices in SL(d, Z) and let D C P(M‘^) be 
adapted. Let i?(®) (respectively Bn{x)) denote the transposed o/T(*) (resp. Anix)). If for all i € A 
we have 


bA 


< 1 . 
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Then for any point x = XqXi ... G A we have 


\\Bn{x)\\D„{x) < 1 - 

In particular, if x = {xqXi .. .Xp-i)°° is periodic the matrix Ap{x) = .. .A^^p-^'> has at 

most one eigenvalue greater than one in absolute value. 

Proof. Let B be the transposed cocycle. The hypothesis is just the case n = 1. Assume that 
this inequality holds for n. By definition An+i{x) = Ai[x)AnlTx) and Dn+i(.x) = Ai(x)Dn(Tx). 
Hence v S DniTx) if and only if Ai{x)v G Dn+i{x). Let us choose v G Dn{Tx), then 

\\Bn+i{x)\\Ai(x)v = \\Bn{Tx)Bi{x)\\A^(x)v < ||H„(ra:)||„ • \\Bi{x)\\ai(x)v 

Now, let X = {xqXi .. .Xp-i)°° be a periodic point. Because Dp_i(x) D Doo{x), the matrix 
Bp{x) satisifies 

\\Bp{x)\\d^(^x) a 1- 

By Lemma|^a basis of the dominant eigenspace of Bp belongs to Dao(x). Consequently, the union of 
the orthogonals of D^oix) contains all eigenspaces corresponding to the non-dominant eigenvalues. 
From the above inequality, we deduce that the absolute value of all eigenvalues different from the 
first one are bounded by 1. □ 


4 Pisot property for Arnoux-Rauzy matrices 

In this section we prove Theorem Let d > 2 be an integer and let ei, 62 , ... be the canonical 
basis of Let e = ei -|- 62 -I- ... -I- and for j = 1,..., fc let /i = e — e^. The domain Dps.d is the 
convex hull of the rays vectors K+/i. 

Let as usual B<-^'> = and (x) = {An{x))*. We claim that we even have a stronger 

property than what is required in Lemma 

Vz = i,...,d, |1 h«||d<i. 

Let us prove this claim. Let v G Dao, then we may write v = pifi + ^ 2/2 -I- ... -f p,dfd for some 
non-negative numbers pLi that satisfy /ii-|-/i 2 -l-...-l-/id = l. We hence have v = e — 

Hy = {z G {z, u) = 0} = {z G ^ Zj = ^ 

Given z G Hy we have B^lz = (zi,..., Zi_i, ^ p^zZj, Zj+i,..., Zd), In other words, B^l acts on Hy 
as a stochastic matrix P{i, v) which is the identity except its Tth row which is {pi, fj. 2 ,..., fid)- In 
particular, < 1 . 

Now for a given finite product A = A^''°'l if one of the letter {1,..., d} is missing 

in the sequence (ioAii ■ ■ ■ Ap-i) then Ae^ = Ci and so the matrix is not primitive. On the other 
hand, if all letters appear it is easy to see that all entries in A are positive. 

Now let X = {xqXi ... Xp-i)°° be a periodic point that contains all letters from A. Because of 
positivity, all the orbits is contained in the interior of D and the dominant eigenvalue is simple. 
Let uo G D be a dominant eigenvector and let = An{x)vQ. Because all (u„) belongs to the 
interior of D the coefficients fii, fi 2 , ..., fid that appear in the stochastic matrices P{xi, vf) are all 
positive. Now, the product P{xp-i,Vp-i)... P{xi,vi)P{xo,vo) is a stochastic matrix with all its 
entries positive. Hence, its second eigenvalue, which is also the second eigenvalue of Ap{x), is less 
than 1 in absolute value. 


5 Pisot property for Brun algorithm (in dimension 3) 

We now turn to the proof of Theorem ,1 Let A = {1,2,3} and A^^\ A^^^ A^^^ be the matrix 
of the Brun algorithm. We let A = and denote by A and B respectively the cocycle and the 
transposed cocycle. We claim that, as in the case of the fully subtractive, we have the stronger 
property that 

VzG {1,2,3}, |1 hW||^„<1. 
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(1 0 0 \ 

We only need to consider the matrix M ^ since the other two are obtained by 

... VO 0 V 

multiplying by a permutation matrix which will not change the L°°-norm. 

Let V = fii{l : 0 : 0 ) + ^2(1 : 1 : 0 ) + ^3(1 : 1 : 1 ) € D for some such that 

fJ-i + 1^2 + fJ'3 = 1 and Hy = {z G + Z2 = /ri^2 — Now, for any z € Hy we have 

= (^1,2:1 + Z 2 ,Z 3 ) = {zi,^iZ2 - z^i). In other words ||S(^)(zi, Z2, 2 :a)l|i < 

II (^ 1 ,2:2,2:3) 11 - 

Now, given a product A = it is easy to see that if 3 does not appear in the 

sequence (* 1 ,^ 2 ) • • • An) then ^63 = 63 and hence the matrix A can not be irreducible. Conversly, 
if 3 appears then A^ is easily seen to be positive. 

Now consider the matrix P built from the begining of the proof. As in the case of the fully 
subtractive algorithm for a primitive product we got that the jii are all positive. Given a product 
A where the matrix appears, the matrix A^ is then such that all rows are such that sum of 
their absolute values is strictly less than one. In other words ||A^|| < 1. 


6 Lyapunov exponents 

Let be a finite or countable set of matrices. Let A, T, A, B denote as before the infinite 

words, the shift map the cocycle and the transposed cocycle. Let also D be adapated to these 
matrices. 

The asymptotic of the cocycle (or the transposed cocycle) are studied through Lyapunov ex¬ 
ponents. Given a T-invariant ergodic probability measure p on A, we associate the real numbers 
li A J 2 — ■■■ — "fd defined by 




7i + 72 + • ■ • + 7fc = lim 

n—yoo 


log II a'® An{x) 


A 


-dp{x). 


In order to be well defined we assume that 


IA 


max (log ||Ai(a;)||,log ||Ai(a:) ^||) dp{x) < co 


( 1 ) 


and we refer to this condition as the log-integrability of the cocycle. If the alphabet A is finite the 
cocycle is automatically log-integrable. If a: is a periodic point of T and p = + 5tx + ... + 

is the sum of Dirac masses distribued along its orbit, then the associated Lyapunov 
exponents are the logarithms of the absolute values of eigenvalues of An{x) where n is the period 
of X. In that sense, Lyapunov exponents generalize eigenvalues. 

Given a measure p for which the cocycle is log-integrable, we say that (A, T, A, p) has Pisot 
spectrum if the associated Lyapunov exponents satisfy 7 ^ > 0 > 72 . This property is related to 
the strong convergence of higher dimensional continued fraction algorithm |La93| . 

Now we restate Lemmain a more dynamical context. 

Lemma 6. Let be a finite or countable set of non-negative matrices in SL(d,Z). Let 

(A,T, A, B) be the associated full shift with its cocycle and its transposed cocycle. Let also D be 
adapated. Assume that 


Vi e A, 




DG) 


< 1 . 


Let p be a T-invariant and ergodic measure on D so that 

• the cocycle A„ is log-integrable, 

• there exists a cylinder [re] such that /idw]) > 0 , A^®") is positive and IjA*^®" ^IIdc™) < 1- 
Then two first Lyapunov exponents of the cocycle An for the measure p satisfies yf > 0 > 72 . 


Proof. Let us first prove that 71 > 0. 

Now, by definition, for /i-almost every x 


7 i = lim 

n—>-oo 


log||A„(a;)|| 


Let m = Iwl be the length of w and consider the position which are multiple of m. For a ^-generic 
X we have by Birkhoff theorem that 

lim fliTAiTTAiM = ,,([„]) 

n—>-oo Tl 
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In other words, given sequence of length n large enough we can find a linear number of disjoint 
occurrences of w (up to a sublinear error). Let kn be the number of these occurrences, then 
necessarily each entry of An{x) is larger than the corresponding one in C*" where C is the matrix 
which contains a 1 in every position. In particular 71 > 0. 

From the existence of w it also follows that for a ^-generic x the cone Doo{x) is reduced to a line 
contained in the interior of KlJ.. We can hence define /i-almost everywhere a function v : A —>■ 
by Doo{x) = K+u(a;) and ||u(a:)|| = 1. We then have the following formulas which holds for /r-almost 
every x 


7 i = lim 

n—foo 


-logP„(a:) ^u(a;)|| 


n 


d^{x) and 


72 = lim 

n—^oo 


log|lS„(a;)||^(2,) 


d^{x). 


It is then easy to derive the estimate for 72 . The map x 1 —)■ v{x) and the dual hyperplanes Hy(^x) 
satisfy the following covariance properties 

Dao{Tx) = A[x)~^Doo{,x) and Ha-i^ = A* H^. 

Hence as in the proof of Lemmawe deduce that for all v S Daoix) 

\\Bm+n{x)\\y(^^) < \\Brn{x)\\y(x)\\Bn{T"^x)\\y{Trr.^) 

In particular, for ^-almost every x £ [re], any n > |?c| we get that \\Bn{x)\\y(^x) < 1- Let S = 
||bW|| 

£)(ra) < 1. Using the same argument as in the estimation of 71 we get that 

, , fcnlog(5 

72 < limini-. 

n—>-oo Ti 

And the above limit is strictly negative. □ 
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